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ABSTRACT 

We  introduce  a random  field  model  with  anisotropic  interfacial  exchange  for  ferromagnet 
/ antiferromagnet  bilayers  and  calculate  average  coupling  energies  and  coupling 
directions  for  different  grain  or  domain  sizes.  The  model  is  shown  to  reconcile 
MalozemofFs  random  field  Ising  model  for  exchange  bias  with  the  spinflop  coupling 
mechanism.  Furthermore,  we  find  that  for  small  and  mid  size  domains,  the  average 
coupling  direction  strongly  varies  between  different  interfacial  configurations.  We  show 
how  this  behavior  can  lead  to  rotating  anisotropies  in  systems  with  very  strong 
magnetocrystalline  anisotropies,  where  the  use  of  superparamagnetic  grains  or  domains 
in  the  antiferrofmagnet  is  not  justified. 


INTRODUCTION 

The  most  frequently  discussed  effects  in  ferromagnetic  (FM)  / antiferromagnet  (AFM) 
bilayers  is  the  shift  in  the  magnetization  curve  that  can  be  observed  when  the  system  has 
been  cooled  in  an  external  ficldf  1 ] . The  magnitude  of  the  shift  is  inversely  proportional  to 
the  thickness  of  the  FM  layer  which  implies  that  the  effect  is  directly  related  to  the 
interfacial  exchange  between  the  FM  and  AFM  and  the  effect  is  thus  referred  to  as 
exchange  anisotropy  or  exchange  bias  [2, 3].  The  temperature  below  which  the  exchange 
bias  is  observed  is  called  blocking  temperature  (TB),  it  is  typically  equal  or  somewhat 
smaller  than  the  Neel  temperature  (Tn)  of  the  AFM.  It  is  commonly  believed  that  the  bias 
sets  in  when  the  AFM  spins  order  and  a sufficient  number  of  domains  in  the  AFM  are 
blocked  such  that  the  AFM  spins  do  not  reverse  when  the  FM  magnetization  is  inverted 
during  the  measurement  of  the  magnetization  curve. 

In  most  cases,  hysteretic  effects  such  as  a strong  increase  in  coercivity  are 
observed  along  with  the  onset  of  exchange  bias[4-6].  Most  models  to  explain  these 
effects  are  based  on  an  idea  of  Fulcommer  and  Charap  [7],  where,  assuming  that  the  FM 
couples  in  some  form  to  the  AFM,  magnetic  structure  of  smaller  grains  or  domains  in  the 
AFM  will  be  dragged  along  irreversibly  with  the  FM  magnetization  and  lead  to 
hysteresis,  while  the  larger  grains,  in  which  the  magnetic  order  is  locked  in,  wifi  lead  to 
the  bias  effect.  Stiles  and  McMichael  [8]  introduced  the  notion  of  a rotating  anisotropy 
which  follows  the  FM  magnetization  when  the  external  field  is  rotated  and  is  attributed  to 
magnetically  “loose”  grains  in  the  AFM.  With  this  model  they  were  able  to  explain  a 
direction  independent  shift  in  the  resonance  field  of  ferromagnetic  resonance 
experiments. 
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Figure  I : Model  of  chree  spins  where  the  F-A  pan  of  the  exchange  energy  is  given  by 

1*1.2  <=1,2 

While  the  concept  of  a superparamagnetic  structure  in  smaller  grains  or  domains 
in  the  AFM  is  reasonable  for  systems  with  weak  magnetocrystalline  anisotropy  such  as 
NiO,  it  is  not  clear  why  is  should  apply  to  an  AFM  like  CoO  where  the  anisotropy  energy 
is  of  the  same  otder  as  the  exchange  energy{9].  In  the  present  contribution  we  show  that 
within  the  random  field  model  a variable  effective  coupling  direction  that  arises  from 
adjustment  of  the  walls  in  smaller  and  mid-sized  AFM  domains  will  lead  to  a rotating 
anisotropy.  In  this  picture  most  AFM  moments  are  frozen  and  dissipative  effects  are  due 
to  irreversible  motion  of  domain  walls. 

RANDOM  FIELD  MODEL  WITH  ANISOTROPIC  EXCHANGE 

In  order  to  understand  how  the  net  coupling  direction  of  the  FM  to  an  AFM  domain  can 
vary  without  rotating  the  AFM  domains,  we  first  have  to  investigate  the  microscopic 
coupling  mechanisms  between  the  AFM  and  the  FM  moment.  For  interfaces  with 
compensated  exchange  between  the  FM  and  the  AFM,  i.e.  where  the  FM  spins  couple 
equally  to  all  AFM  sublattices[10],  the  following  two  models  are  typically  discussed  in 
theliteraturefll,  12]: 

1.  The  random  field  model  [13-15],  in  which  the  AFM  is  assumed  to  break  up  into 
domains.  The  finite  area  of  these  domains  will  result  in  a net  interfacial  exchange 
which  scales  as  (i/Jn)  if  N is  the  number  of  AFM-FM  pairs.  In  this  model  the 
FM  and  AFM  spins  are  assumed  to  be  collinear  which  leads  to  a change  in 
exchange  energy  A£  = 27  /-Jn  when  the  FM  spins  are  reversed.  The  random 
field  model  is  commonly  used  to  explain  the  net  unidirectional  coupling  between 
the  FM  and  individual  AFM  domains  (or  grains). 

2.  The  spin-flop  coupling  model  [ 1 6, 1 7],  which  assumes  Heisenberg  spins  and 
where  the  frustrated  exchange  across  the  interface  leads  to  a slight  canting  to  the 
AFM  moments  and  a net  coupling  of  the  FM  that  is  perpendicular  to  the  AFM 
easy  axis.  In  contrast  to  the  random  field  Ising  model  tire  coupling  is  uniaxial  and 
does  not  lead  to  exchange  bias  [18, 19].  Furthermore,  the  coupling  strength  is 
independent  of  the  AFM  domain  / grain  size. 

The  coDinearity  required  by  the  random  field  Ising  model  is  only  justified  for  small 
domains  in  the  AFM.  For  most  realistic  domain/grain  sizes,  the  spin-flop  mechanisms 
will  dominate  [8]  leaving  the  origin  of  the  unidirectional  coupling  as  a puzzle.  We  will 
now  show  that  spin-flop  coupling  can  have  a unidirectional  component  when  the 
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anisotropic  Dzyaloshinsky-Moriya  (DM)  term  is  added  to  the  exchange,  i.e.  the  exchange 
energy  between  two  moments  on  sites  i and  j has  the  form 

EU  =-J,js,si+D:j{sixIi).  (1) 

In  cases  where  the  crystal  structure  has  a center  of  inversion  at  the  point  bisecting  the  line 
between  i and  j,  Z>  can  be  shown  to  vanish  [20, 21].  While  this  symmetry  condition 

holds  in  most  bulk  system,  it  is  never  satisfied  at  the  AFM/FM  interface[22]  and  the  DM 
term  thus  has  to  be  included  to  properly  describe  interfacial  exchange.  The  consequence 
of  the  DM  term  is  that  the  difference  in  exchange  energy  upon  reversal  of  the  FM  spins  is 
different  from  zero  in  the  spin-flop  coupling  state.  This  is  best  illustrate  by  considering 
three  spins  (see  Figure  1),  SF,  SAJ , and  SA2 , in  the  spin-flop  coupling  states.  The  energy 
difference  when  SF  is  reversed  (Figures  la  and  lb)  is  AE^  = 2(d,  . + Dlr ) . 

In  order  to  estimate  the  magnitude  of  the  unidirectional  term  in  the  case  of  a 
FM/AFM  bilayer,  we  perform  a calculation  using  the  random  field  model  with 
anisotropic  interfacial  exchange.  We  consider  a bilayer  with  bcc  type  lattice  and  a (1 10) 
interface  with  nearest  neighbor  exchange  / in  the  FM  and  -J  in  the  AFM.  In  the  AFM  we 
further  assume  uniaxial  anisotropy  along  the  (001)  direction  with  KU~  1 . For  the 
interfacial  exchange  we  use  Eq.  (1)  with  Jtj  = Jr  and 

Dy = Dr' (sin dcostp, sin tfsin^, cost?),  where  re  [-1,1],  r'e  [-1,1],  tie  ]0,ff) , and 
<p  6 [0,2/r)  are  uniform  random  variables.  These  random  fluctuations  are  usually  justified 
with  the  large  variation  of  the  interfacial  bond  lengths  that  results  from  the  large  lattice 
mismatch  between  the  FM  and  AFM  in  systems  like  CoO/Permally  or  FeF2/Fe.  The 
particular  system  we  consider  here  consists  of  5 layers  in  the  AFM  and  10  layers  in  the 
FM.  We  use  periodic  boundary  conditions  parallel  to  the  interface  with  variable  unit  cell 
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Figure  2:  Energy  change  due  to  reversal  of  the  FM  spins  for  AFM-FM  bilayers  (see  texl  for  details)  for 
Heisenberg  spins  with  random  anisotropic  interfacial  exchange  and  different  sizes  of  the  computational 
ceil  (symbols)  compared  to  expected  result  from  the  random  field  Ising  model  (solid  line) 
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size  consisting  of  L by  L primitive  cells  ( L = 2, 4, 8, 16,  and  32).  The  effect  of  the 
magnetostatic  interaction  is  treated  as  an  effective  hard  axis  in  the  FM  perpendicular  to 
the  interface.  We  seek  to  calculate  AFM/FM  coupling  energies  as  a function  of  domain 
size,  where  we  equate  the  unit  cell  size  of  our  model  to  the  domains  size.  For  this  we 
have  to  suppress  the  formation  of  further  domains  within  the  unit  cell.  We  have  therefore 
prepared  an  ordered  AFM  spin  configuration  and  randomized  FM  spins  as  an  initial 
condition  and  subsequently  minimized  the  energy  using  the  Landau  Lifshitz  equations 
(see  [ 1 9]  for  details).  When  the  energy  minimum  is  reached,  the  FM  spins  are  inverted 
and  the  energy  is  again  minimized.  For  all  but  the  smallest  unit  cells  this  procedure 
always  yields  two  distinct  states.  Their  energy  differences,  averaged  over  40  random 
configurations,  are  plotted  in  Figure  2 and  compared  to  the  result  of  the  random  field 
Ising  model.  While  the  unidirectional  anisotropy  due  to  the  DM  term  is  systematically 
smaller  than  that  of  die  random  field  Ising  model,  it  has  nevertheless  comparable  order  or 
magnitude  provided  D/J  is  not  too  small.  Moria  [21]  estimated  DU  ~ 0.1  in  bulk  AFM 
but  values  that  are  as  large  as  0.5  have  been  reported  for  systems  with  lower  symmetry 
[23]. 

EFFECTIVE  COUPLING  DIRECTION  AND  ROTATING  ANISOTROPIES 

fa)  the  previous  section  we  have  determined  energy  differences  averaged  over  many 
random  configurations  of  the  interfacial  exchange  interactions.  The  direction  of  the  FM 
magnetization  with  regard  to  the  AFM  easy  axis  for  each  of  these  configurations  is 
plotted  in  Figure  3.  By  choosing  appropriate  initial  conditions  for  the  minimization 
procedure  discussed  above,  the  system  can  always  be  prepared  in  a way  that  the  global 
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figure  3:  Average  direction  of  FM  spins  for  every  random  configuration.  The  AFM  easy  axis,  (001 ). 
is  taken  to  be  parallel  to  the  y-axis  of  the  plot  In  the  case  ZM)  and  i=2  (4),  80%  (50%)  of  die 
configurations  are  collinear  with  the  AFM  easy  axis. 


r-  r . 1 t i " t > r 1 ■ r 1 t 

■ A.A. 

i mmmmk 

I J 1 1 1 - X I I ■ t 1 - 


34 


energy  minimum  has  the  FM  moments  pointing  at  angle  between  0 and  90  degrees  with 
regard  to  the  AFM  easy  axis  (this  choice  of  initial  condition  corresponds  to  the  field 
cooling  procedure).  As  expected,  the  coupling  is  always  collinear  for  small  cells  and 
perpendicular  when  the  cells  are  large.  When  the  anisotropic  exchange  is  non  zero,  the 
coupling  direction  strongly  varies  between  different  configurations  of  the  smaller  and 
mid  size  cell  but  is  always  perpendicular  to  the  AFM  easy  axis  when  the  cells  are  large. 

As  was  done  in  the  previous  section,  we  consider  the  cell  size  in  the  model  to 
correspond  to  domain  sizes  in  the  film  and  further  recognize  that  different  domains 
correspond  to  different  random  configurations  of  the  interfacial  exchange  in  the 
computational  cell.  To  simplify  the  argument,  let  us  assume  that  the  field  was  applied 
along  the  (1-10)  direction  during  the  field  cooling  procedures,  i.e.  perpendicular  to  the 
AFM  easy  axis.  The  domain  pattern  in  die  AFM  will  arrange  itself  such  that  the  coupling 
direction  of  most  domains  is  parallel  to  the  cooling  field.  When  the  direction  of  the 
applied  field  and  the  FM  magnetization  is  subsequently  rotated  to  a different  direction, 
the  AFM  domain  pattern  will  adjust  by  moving  domain  walls  to  minimize  the  energy. 
While  it  remains  to  be  shown  explicitly,  the  displacement  of  such  domain  walls  in  real 
samples,  with  defects  at  which  these  walls  are  pinned,  is  likely  an  irreversible  process. 
According  to  the  results  plotted  in  Figure  3,  the  coupling  direction  of  the  largo'  domains 
will  not  be  affected  by  the  rearrangement  of  domain  walls.  Together  with  the  energy 
differences  plotted  in  Figure  2,  this  will  lead  to  a unidirectional  anisotropy  that  has 
similar  order  of  magnitude  as  that  found  with  the  random  field  Ising  model.  Since  the 
coupling  direction  of  the  smaller  domains  depends  strongly  on  the  random  configuration, 
changing  the  domain  pattern  will  lead  to  new  configuration  with  adjusted  coupling 
directions  that  follows  the  FM  magnetization.  This  will  give  rise  to  a rotating  component 
of  the  anisotropy  which  adjusts  with  direction  of  the  applied  field,  hi  this  model  the 
dissipation  would  take  place  in  the  irreversible  displacement  of  domain  walls  in  the  AFM 
and  thus  does  not  require  the  superparamagnetic  rotation  of  all  the  moments  in  smaller 
grains. 

ACKNOWLEDGMENTS 

This  research  was  sponsored  by  the  Laboratory  Technology  Research  Program  of  the 
Work  supported  by  Department  of  Energy,  Office  of  Science  through  BES-DMSE,  and 
OASCR-M1CS  as  well  as  the  Laboratory  Technology  Research  Program  of  OASCR- 
TRD  under  Contract  No.  DE-AC05-00OR22725  with  UT-Battelie  LLC. 

REFERECES 

1 . Meiklejohn,  W.H.  and  C.P.  Bean,  A new  magnetic  anisotropy.  Physical  Review 
B,  1956. 102:  p.  1413. 

2.  Nogues,  J.  and  LK.  Schuller,  Exchange  bias.  Journal  of  Magnetism  and  Magnetic 
Materials,  1999. 192(2):  p.  203-232. 

3.  Berkowitz,  A.E.  and  K.  Takano,  Exchange  anisotropy  - a review.  Journal  of 
Magnetism  and  Magnetic  Materials,  1999. 200(1-3):  p.  552-570. 

4.  Stiles,  M.D.  and  R.D.  McMkhael,  Coercivity  in  exchange-bias  bilayers.  Physical 
Review  B,  2001 . 6305(6):  p.  art.  no.-064405. 

5.  Li,  Z.  and  S.  Zhang,  Coercive  mechanisms  in  ferromagntic-antiferromagnetic 
bilayers.  Physical  Review  B,  2000. 61(22):  p.  R14897-R14900. 


35 


6.  Leighton,  C.,  et  al„  Thickness-dependent  coercive  mechanisms  in  exchange- 
biased  bilayers.  Physical  Review  B,  2002. 65(6):  p.  art.  no. -064403. 

7.  Fulcomer,  E.  and  SJL  Charap,  Thermal  Fluctuation  Aftereffect  Model  for  Some 
Systems  with  Ferromagnetic-Antiferromagnetic  Coupling.  Journal  of  Applied 
Physics,  1972. 43(10):  p.  4190-&. 

8.  Stiles,  M.D.  and  R.D.  McMichael,  Model  for  exchange  bias  in  potycrystalline 
ferromagnet-antiferromagnet  bilayers.  Physical  Review  B,  1999. 59(5):  p.  3722- 
3733. 

9.  Wagner,  V.  and  D.  Hermann-Ronzaud,  Neutron  Inelastic  Scattering  1977.  IAEA, 
Vienna  1977, 1977.  Part  H:  p.  135-143. 

10.  Here  we  will  discuss  uniaxial  AFMs  that  consist  of  two  ferromagnetic  sublattices 
with  anti-aligned  moments  and  we  assume  nearest  neighbor  exchange  between 
AFM  and  FM  spins  at  the  interface  with  exchange  parameter  J 

11.  Ki  wi,  M.,  Exchange  bias  theory.  Journal  of  Magnetism  and  Magnetic  Materials, 
2001. 234(3):  p.  584-595. 

12.  Stamps,  RJL,  Mechanisms  for  exchange  bias.  Journal  of  Physics  D- Applied 
Physics,  2000. 33(23):  p.  R247-R268. 

1 3.  Malozemoff,  A.P.,  Random-Field  Model  of  Exchange-Anisotropy  at  Rough 
Ferromagnetic-Antiferromagnetic  Interfaces.  Physical  Review  B,  1987. 35(7):  p. 
3679-3682. 

14.  Malozemoff,  A.P.,  Mechanisms  of  Exchange-Anisotropy.  Journal  of  Applied 
Physics,  1988. 63(8):  p.  3874-3879. 

1 5.  Nowak,  U.,  et  al..  Domain  state  model  for  exchange  bias.  I.  Theory.  Physical 
Review  B,  2002. 66(1):  p.  art.  no.-014430. 

16.  Koon,  N.C.,  Calculations  of  exchange  bias  in  thin  films  with 
ferromagnetic! antiferromagnetic  interfaces.  Physical  Review  Letters,  1997. 
78(25):  p.4865-4868. 

1 7.  Hinchey,  L.L  and  D.L  Mills,  Magnetic-Properties  of  Ferromagnet- 
Antiferromagnet  Superlattice  Structures  with  Mixed  Spin  Antiferromagnetic 
Sheets.  Physical  Review  B,  1986. 34(3):  p.  1689-1699. 

18.  Schuhhess,  T.C.  and  W.H.  Butler,  Consequences  of  spin-flop  coupling  in 
exchange  biased  films.  Physical  Review  Letters,  1998. 81(20):  p.  4516-4519. 

19.  Schulthess,  T.C.  and  W.H.  Butler,  Coupling  mechanisms  in  exchange  biased  films 
(invited).  Journal  of  Applied  Physics,  1999. 85(8):  p.  5510-5515. 

20.  Dzyaloshinsky,  L,  A thermodynamic  theory  of  "weak"  ferromagnetism  of 
antiferromagnetics.  J.  Phys.  Chem.  Solids,  1958. 4:  p.  241-255. 

21.  Moriya,  T.,  Anisotropic  Superexchange  Interaction  and  Weak  Ferromagnetism. 
Physical  Review,  1958. 120(1):  p.  91-98. 

22.  Crepieux,  A.  and  C.  Lacroix.  Dtyaloshinsky-Moriya  interactions  induced  by 
symmetry  breaking  at  a surface  Journal  of  Magnetism  and  Magnetic  Materials. 
1998.182(3):  p.341-349. 

23.  Katsnelson,  M.L,  V.V.  Dobrovitski,  and  B.N.  Harmon,  Many -spin  mode l and  the 
spin  Hamiltonian  ofMn-12  clusters.  Journal  of  Applied  Physics,  1999. 85(8):  p. 
4533-4535. 


36 


